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Abstract. Here we present the results of panoramic and long-slit observations of eight ULX 
nebular counterparts held with the 6m SAO telescope. In two ULXNe we detected for the first 
time signatures of high excitation ([0 iii]A5007 / H/3 > 5). Two of the ULXs were identified with 
young (T ~ 5 — 10 Myr) massive star clusters. Four of the eight ULX Nebulae (ULXNe) show 
bright high-excitation lines. This requires existence of luminous (~ 10 38 -M0 40 ergs s" 1 ) UV/EUV 
sources coinciding with the X-ray sources. Other 4 ULXNe require shock excitation of the gas 
with shock velocities of 20-100fcm s _1 . However, all the studied ULXN spectra show signatures of 
shock excitation, but even those ULXNe where the shocks are prevailing show presence of a hard 
ionizing source with the luminosity at least ~ 10 38 ergs s _1 . Most likely shock waves, X-ray and 
EUV ionization act simultaneously in all the ULXNe, but they may be roughly separated in two 
groups, shock-dominated and photoionization-dominated ULXNe. The ULXs have to produce 
strong winds and/or jets powering their nebulae with ~ 10 39 ergs s _1 . Both the wind/jet activity 
and the EUV source needed are consistent with the suggestion that ULXs are high-mass X-ray 
binaries with the supercritical accretion disks of the SS433 type. 
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1. Introduction 

1.1. Ultraluminous X-ray Sources 

A point-like X-ray source is considered an Ultralu- 
minous X-ray Source (ULX) if its luminosity exceeds 
10 39 erg s^ 1 and it does not coincide with the galac- 
tic nucleus (to exclude active galactic nuclei). There 
are more than 150 ULXs known at the present time 
(Swartz, 2004) but very little is clear yet about the 
physical nature of these objects. 

Luminosity of an accreting source is limited by the 
Eddington value, that for a stellar mass (~ 10M Q ) 
black hole is about 10 39 ergs s _1 . The limit, how- 
ever, can be violated by a logarythmic factor if ac- 
cretion is supercritical (Shakura & Sunyaev, 1973; 
Abramowicz et al., 1980). ULXs most likely are ac- 
creting black holes (their compactness follows from 
relatively fast variability of some of them (Krauss et 
al., 2005; Soria et al., 2006)), and their apparent lu- 
minosities are well above the Eddington limit for a 
stellar mass (~ 10M Q ) black hole. Different models 
were proposed to explain the ULX phenomenon, bas- 
ing on three main ideas: 

- the black hole mass may be larger, ~ 100 — 
10000 solar masses, - Intermediate Mass Black Holes, 
IMBHs (Madau & Rees, 2001; Colbert & Miller, 2005) 

- accretion may be supercritical but still radia- 



tively efficient (Begelman, 2002) 

- observed source may be intrinsically 
anisotropic (Fabrika & Mescheryakov, 2001; King et 
al., 2001) 

IMBHs are supposed to accrete matter from mas- 
sive donor stars, accretion from interstellar medium 
is unlikely to provide enough material. Possible obser- 
vational appearences of IMBH+donor close binaries 
were considered by Hopman et al. (2004) and Cop- 
perwheat et al. (2005). IMBHs are considered either 
compact remnants of Population III stars (Madau & 
Rees, 2001) or products of stellar collisions in stellar 
cluster or protocluster cores (Soria, 2006). 

Supercritical accretion disks (SCADs) are the 
most appreciated alternative to the scenarios involv- 
ing IMBHs. The best candidates are massive black- 
hole binaries like SS433, the only example of a super- 
critical accretor in the Galaxy (Fabrika, 2004). 

1.2. Optical Counterparts 

For most of the well-studied ULXs optical counter- 
parts were found. In some cases optical observations 
reveal stellar counterparts - usually OB supergiants 
(Liu et al., 2001; Terashima et al., 2006). Many of 
ULXs are situated inside star-forming regions and 
therefore they are heavily absorbed (Terashima et al., 
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2006). 

Large number of ULXs are located inside neb- 
ulae (ULX nebulae, ULXNe), sometimes superbub- 
bles more than hundred parsecs in size (Pakull & 
Mirioni, 2003; Pakull et a!., 2006), usually having 
spectra typical for shock-powered nebulae (bright 
[S n]A6717,6731, [N n]A6748,6583, [O i]A6300,6364 
and [N i]A5200 lines). In some cases the appearence 
of a ULXN may be quite different, like the high- 
excitation photoionized H n region of HoII X-l 
(Lehmann et al., 2005) or a compact very bright shell- 
like nebula with anomalous excitation conditions (the 
nebula MF16, (Abolmasov et al.2006; Abolmasov et 
al.2007a)). 

One of the interesting properties of some ULXNe 
is He nA4686 emission line. Due to this rea- 
son, ULXNe were considered X-ray Ionized Nebu- 
lae (XINe, Pakull & Mirioni, 2003). This is often 
considered an argument for truely large X-ray lumi- 
nosities of ULXs (Kaaret et al., 2004). However, ex- 
treme ultraviolet radiation in the wavelength range 
A ~ 100 — 200 A has larger ability for producing 
He IIA4686 recombination emission line than X-rays. 
Therefore nebulae like MF16 and that around HoII X- 
1 may be not X-ray Ionized Nebulae. For more details 
see discussion in Abolmasov et al. (2007a). In some 
cases like HoIX X-l and NGC1313 X-l (Grise et al., 
2006; Pakull et al., 2006 ) the Hell emission is broad- 
ened by ~ 1000/cto s" 1 and can be produced either 
in the accretion disk or in a Wolf-Rayet donor atmo- 
sphere. 

Many of ULXs are located close to young stellar 
clusters, often at tens to hundreds parsecs from the 
clusters or groups of massive stars (Zezas et al., 2002). 
There are numerous indications for the ULXs are re- 
lated to young (5-10 Myr) stellar population (Soria et 
al., 2005; Zezas et al., 2002; Abolmasov et al., 2007b). 
For distant ULXs (D > 15 Mpc) bright star clusters 
with surrounding nebulae become the only associa- 
tion available for optical studies. 

In the following section [2] we describe our obser- 
vations of ULX counterparts with the 6m BTA tele- 
scope. In section [3] results of the spectral analysis and 
modelling of the optical counterparts are presented. 
In section @] we discuss the inference on the nature of 
ULXs which may follow from the observations of the 
ULXNe. 

2. Observations and Data Reduction 

We have focused on targets with known optical coun- 
terparts, situated in nearby galaxies (D < lOMpc). 
The only exception is NGC7331 X-l, but in that case 
a bright host cluster of V ~ 20 m is present. All the 
host galaxies are of late Hubble types. In table [2] we 
list the ULX counterparts observed. We used accu- 
rate Chandra X-ray coordinates (±0.3") according to 



Swartz et al.(2004), Roberts et al.(2003) in the case 
of IC342 X-l, Terashima & Wilson (2004) in the case 
of M51 X-7 and Mescheryakov (2004) for HoIX X-l. 

All the data were obtained with two spec- 
trographs: Multi-Pupil Fiber Spectrograph MPFS 
(Afanasiev et al, 2001) and SCORPIO (Afanasiev 

6 Moiseev, 2005) focal reducer in long-slit mode. 3D 
data have the advantage of providing unbiased flux 
calibrations. However, SCORPIO has higher quan- 
tum efficiency by a factor of up to ~ 6. 

MPFS was used with the grating # 4 
(600 lines/mm), SCORPIO was used with grisms 
VPHG550G and VPHG1200G. Spectral ranges and 
spectral resolutions are listed in table [2] Seeing con- 
ditions during observations were from 1 to 2" . For 
the long-slit observations we estimate the light losses 
on the slit, consistent with the actual seeing value. 
For the sources listed in table [2] we have obtained 
either long-slit or panoramic spectra appropriate for 
analysis. 

Data reduction system was written in IDL6.0 en- 
vironment, using procedures written by V. Afanasiev, 
A. Moiseev and P. Abolmasov. For long-slit observa- 
tions the reduction system contains all the standard 
procedures. For 3D data the standard procedures for 
panoramic data reduction were used (i. e. those stan- 
dard for long-slit data plus extraction of fiber spec- 
tra and fiber sensitivity calibration) . For MPFS data 
atmospheric dispersion correction was made with ac- 
curacy ~ 0".l. 

In all the targets emission lines from the sur- 
rounding nebulae were detected. Contribution from 
the stellar continua (host cluster or association) has 
also been detected in some cases. Emission line pa- 
rameters were calculated using gauss analysis. Fitting 
by two gaussians was used for [SII]AA6717,6731 dou- 
blet, H 7 +[OIII]A4363, Hell A4686+FeIIIA4658 and 
Fell+Felll A5262,5270 blends, triple gaussian was 
used to deblend Ha with [Nil] A6548, 6583 doublet (for 
the components of the latter both fixed wavelength 
difference and flux ratio F(A6583/A6548) = 3 were 
adopted). Exact wavelength values with an accuracy 
of ~ Q.1A were taken from Coluzzi (1996). 

Among our targets there are two large-scale bub- 
ble nebulae (HoIX X-l and IC342 X-l) and four 
high-excitation nebulae (HoII X-l, NGC6946 ULX-1, 
M101 P098 and NGC5204 X-l) In two cases (M51 X- 

7 and NGC7331 X-l) comparatively high S/N cluster 
spectra with nebular contribution where obtained. 

3. Results and Interpretation 

Emission line integral fluxes and luminosities are 
given in table [2l In the last row we estimate the to- 
tal power of the nebula suggesting it either a shock- 
powered (see section l3Jj) or photoionized HII region 
(see section I3.2|) . 



Table 1: Observed ULXs. The columns present the target IDs, ULX X-ray coordinates (JD 2000), distances D to the host galaxies in Mpc, Galactic absorption 
Ay in the source direction according to Schlegel (1998), angular distances R of targets towards the center of their galaxies in arcminutes, oxygen abundances 
according to Pilyugin et al. (2003) (the central abundance, but interpolated abundances to the target angular distance from the center are given in brackets), 
comments on the optical counterparts, exposure time in seconds, spectral range in A, spectral resolution in A, in the MPFS oservations and the same 
parameters in the SCORPIO observations. 
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8.85 (8.25) 
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10 
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(Blair et al., 2001) 
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1998) 
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a Saha et al.(2002) 

b Karachentsev et al.(2002) 

c Tully (1988) 

d Tully et al.(1992) 

e Feldmeier et al.(1997) 

f Stetson et al.(1998) 

9 Hughes et al.(1998) 

h Makarova et al.(2002) 



Table 2: Observed HP fluxes, emission-line fluxes in HP units and unreddened luminosities of the ULXNe. For NGC6946 X-l, HoII X-l, M101 P098 and 
NGC5204 X-l the luminosities were unreddened using Ha/Hp criterion and Cardelli et al.(1998) reddening curves. For NGC7331 X-l the best-fit interstellar 
absorption (Ay = 1™A) was derived from the stellar population spectrum fit. In other ULXNe the galactic Ay was used (compare with tabled In NGC7331 X-l 
the luminosities in residual emissions are given in brackets (see text for details). The total power estimates for HoII X-l, NGC6946 X-l, M101 P098 and 
NGC5204 X-l are hydrogen-ionizing luminosities derived in assumption that HP is a recombination line (see section \3.2]) . Remaining nebulae are treated as 
shock-powered, their total powers were estimated using equation (0). 
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3.1. Bubble nebulae 

We classify HoIX X-l MH9/10 nebula and IC342 X- 
1 nebula as purely shock-powered nebulae, with the 
shock velocities in the range 20 — lOOfcm s _1 . Both es- 
timates follow from the kinematical data, however, in 
this paper we concentrate on emission line intensities 
and diagnostics, leaving kinematical studies to fur- 
ther papers. If one suggests photoionisation for these 
two ULXNe, a source with L > 10 42 ergs s _1 is re- 
quired. Their positions on ionization diagrams (see 
Fig. 2]) are consistent with shocks of moderate power 
(Vs ~ lOOfcms -1 ). The total power of a nebula and 
the mean shock velocity (in the case when the physi- 
cal size is known) can be estimated using its total H/3 
luminosity. According to Dopita & Sutherland (1996), 
total energy flux in a shock wave: 

F tot = 2.28 x 10~ 3 U 2 3 n erg cm^s- 1 (1) 

where Vi is the shock velocity in 100 fcms -1 units, 
tiq is the pre-shock hydrogen density in cto~ 3 . For 
the energy flux in H/3: 

^Hfi — ^Hf3, shock + ^H0^precursor 

(7.44 x 10- 6 V 2 2 - 41 + (2) 
+9.86 x 10~ 6 V^- 28 ) n erg cm- 2 s- 1 , 

Fh/3, shock represents the emission of the cooling ma- 
terial behind the shock, Fjjp, precursor is the contri- 
bution of the precursor. The Fn/3,shock term is ap- 
propriate for any steady one-dimensional shock but 
the precursor term goes to zero at Vs < 150km s^ 1 . 
Combining (fl]) and (J3J) one expresses the H/3 luminos- 
ity as a function of the total shock power and shock 
velocity: 

L Hl3 = (3.26 x lO- 3 ^°- 59 + 

+4.32 x 10- 3 V 2 -°- 72 Q(V 2 - 1.5)) L tot [6) 

where Q(x) = 1 if x > and Q(x) — otherwise. 
From this equation the total power can be roughly 
estimated as L tot ~ (100 + 3OO)Lh0- 

In SCORPIO spectra of MH9/10 we detect 
roughly a point-like He nA4686 emission (figure [1} 
close to the blue knot (central star cluster of the 
bubble) and the X-ray source (Grise et al., 2006). 
In figure [2] we present ID maps of the X-ray source 
vicinity demonstrating that the He II line is emit- 
ted at the edge of the central cluster. We estimate 
the luminosity of the He nA4686 source as ~ (2 + 
4)10 36 ergs s _1 (the obtained value was multiplied by 
a factor 2 — 4 estimating the losses of the slit). This 
is consistent with the results obtained by Grise et 
al. (2006), who also detect the He n emission. 

In HoIX X-l nebula this line may have stellar 
origin. Possible sources of bright He IIA4686 lines 
are Wolf-Rayet stars (Conti et al., 1983) and ac- 
cretion disks. For WN stars the equivalent width 
can reach 400A. V-band luminosity is about 10 37 + 
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Figure 2: ID maps for MH9/10 nebula (HoIX X-l) 
derived along the slit in continuum (solid line; inte- 
grated over the wavelength range 4000 — 7000^4, emis- 
sion lines excluded by median smoothing), Ha (dashed 
line) and He 11X4686 (asterisks with error bars). 

10 38 ergs s _1 for bright WN stars, so one WN star can 
be enough to explain the He liA4686emission. 

Bright emission lines in the IC342 X-l nebula 
spectrum argue for shock excitation. However, its in- 
tegral spectrum (figure [3]) containes relatively faint 
high-excitation lines like He 11 and [Fe ill] emitted in 
the inner parts of the nebula. High-excitation spec- 
trum resembles that of MF16 (see next section). 

3.2. NGC6946 X-l 

NGC6946 X-l is known to coincide with a peculiar 
nebula MF16, for a long time considered a lumi- 
nous Supernova Remnant (SNR). It is the most com- 
pact (20pcx34pc) and one of the brightest among the 
ULXNe. Its integral spectrum is characterised by the 
highest [O m]A5007/H/3 and [N n]A6583/Ha ratios. 

The velocity gradients and non-gaussian structure 
of the emission lines in our spectra suggest expan- 
sion velocity not greater than 200fcms~ 1 . Dunne et 
al.(2000) have resolved Ha, [N 11] and [O in] lines in a 
narrow (FWHM^20-40fcm s" 1 ) component from un- 
shocked material and a broad (FWHM~250fcm s _1 ) 
from shocked material. Basing on the MPFS and 
SCORPIO spectra we performed a thoroughful anal- 
ysis of the ionization and excitation sources of the 
nebula (Abolmasov et al., 2006; Abolmasov et al., 
2007a), leading to a conclusion that most of the 
power of the nebula comes not from shock waves 
but from photoionizing ultraviolet source even more 
luminous than the central X-ray source. CLOUDY 
modelling gives the parameters of the central source 
in the suggestion that the ionizing continuum in the 
EUV/UV is a black body: T = io 515±015 X, L BB = 
(7.5 ± 0.5) x 10 39 ergs s" 1 . For the EUV luminosity 
one can make rough, but universal estimates (the 
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Figure 1: MH9/10 (the nebula associated with HoIX X-l) Ha image (taken with SCORPIO) with the 1" slit. 
The X-ray source coordinates are marked by a white circle. To the right the spectra of the two parts of the 
nebula known as MH9 (bottom-right panel) and MH10 (top-right panel). In the upper spectrum He 11X4686 
emission was detected. 



lower limits), based on the ionizing quanta numbers. 
If He nA4686 and H/3 are purely recombination lines 
one can apply the relations by Osterbrock (1974): 



L(X < 228A) > 



4Ry a B (He + + ) 
£(A4686) 

v ' ^Hell A4686 



x L H 



el IX 4686 



may be significant. Discussion about the nature of the 
EUV photoionizing source follows in section [4T21 

High signal-to-noise spectrum of MF16 obtained 
with SCORPIO shows numerous high-excitation lines 
like [Fe in] and Ar IV. Integral spectrum is shown in 
( 4 ) figure El 



100L HeIIX 



II, M> 



L(X < 9121) > 



lRy a B (H+) 



E(Hp) 



xL 



Hp 



65L Hl3 , (5) 



Here as values are total recombination rates for 
the Case B (optically thick to ionizing quanta) 
for hydrogen and ionized helium, correspondingly. 
a e ff (line) is the effective recombination rate with 
the particluar emission line quanta production. 



«s(- ffe++ )/«He / //A4686 



and a.B(H + )/a e I jp ratios de- 



pend weakly on the plasma parameters. The former 
changes by about factor of two in the temperature 
range (0.3 -r 10) x 10 4 X, the latter by about 20%. 

In table [2] we give the hydrogen-ionizing luminosi- 
ties defined this way as estimates for the power of 
photoionized nebulae. H/3 line can overestimate the 
photoionizing flux because contribution from shocks 



3.3. High-Excitation ULXNe 

In the vicinity of two of the observed sources 
(M101 P098 and NGC5204 X-l) we detect regions 
with high [O m]A5007/H/3 ratio (> 5). They can be 
seen as bright spots on the [O in]A5007 flux maps 
shown in figure [6] The [O m]A5007/H/3 ratios are 
similar to that of MF16, but the integral luminosi- 
ties are much less - an order of magnitude less for 
NGC5204 X-l and even more for M101 P098. 

M101 P098 nebula is well resolved as an elongated 
region 1" x 4" , bright in [O m]A5007,4959 lines (see 
figure^)- For the distance 7.2 Mpc this corresponds 
to the physical size about 35pc x 140pc. 

Bright [O ill] emission lines argue for either a pho- 
toionizing source capable to provide ionization pa- 
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Figure 3: Spectrum of the IC342 X-l nebula (SCORPIO). The He 11X4686 emission originates roughly in the 
center of the ULXN. 



rametciQ > 1(T 3 at about 50 pc or shock waves with 
Vs > 300fcm s _1 propagating at similar spatial scales 
(Evans, 1999; Dopita & Sutherland, 1996). Photoion- 
izing source must have the hydrogen-ionizing lumi- 
nosity > 2 x 10 38 ergs s -1 , so M101 P098 can be 
considered an analogue of MF16 in a more rarefied 
ISM. From the other hand, fast shocks with Vs ~ 
300 — 500fcm s _1 require practically identical power 
(1.5 — 2) x 10 38 ergs s _1 to provide the observed H/3 
luminosity. The example of MF16 suggests that these 
high-excitation nebulae are more likely to be pho- 
toionizcd HII regions. However, both photoionization 
and shock excitation can act in powering these nebu- 
lae. 

HoII X-l is another example of a high-excitation 
ULX nebula. [O m]A5007/H/3 - 3 is low if com- 
pared with the M101 P098 and NGC5204 X-l 
spectra, but HoII has lower oxygen abundance, 
so the physical conditions may be quite similar. 
He nA4686 emission line is very bright relative to H/? 



1 We define ionization parameter as (for example, in Evans et 
al.(1999)): 

U = — / — dE, 

cn Jl3.6eV E 



/13.6eV 

where n is the gas number density 



(F(HeIIXA686)/F(H(3) ~ 0.2), but the total lumi- 
nosity in this line is an order of magnitude lower than 
that in MF16. Our results on this object were pub- 
lished in Lchmann et al.(2005). Long-slit spectrum 
(obtained with UAGS spectrograph) of the nebula is 
shown in figure [7) 

3.4. Star Clusters 

Among the observed ULXs two objects coincide 
with isolated bright star clusters: NGC7331 X-l and 
M51 X-7. Among the other ULXs whose X-ray boxes 
contain stellar clusters, there are HoIX X-l (the clus- 
ter is very faint, V ~ 20™5 and not massive) and 
NGC5204 X-l (the cluster is a part of a more com- 
plex star- forming region). 

A nearly point-like continual source is seen in 
ground-based images of the NGC7331 X-l region. 
At the distance of 15 Mpc its size is ~ 50 pc. In 
the HST images the source consists of young clus- 
ter and some fainter sources around (Abolmasov et 
al., 2007b). The cluster coincides with the HII re- 
gion P98 (Petit, 1998). Our study of the source envi- 
ronment based on HST and SCORPIO observations 
is published in a separate paper (Abolmasov et al., 
2007b). We fitted the integral spectrum of the stel- 
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lar population with StarBurst99 models (Leitherer 
et al., 1999; Vazquez & Leitherer, 2005) with solar 
and 0.4 of solar metallicities, and derived the best- 
fit age T — 4.25 ± 0.5 Myr and the cluster mass 
M ~ 4OOOOM . The integral spectrum together with 
StarBurst99 best-fit stellar population and CLOUDY 
nebular spectra is shown in figure [S] WC features are 
bright in the spectrum, pointing to the presence of 



several WC stars. Emission line spectrum can not 
be totally explained by photoionization by the cen- 
tral cluster. Low-excitation shock-ionized lines like 
[S n]A6717,6731 and [N il]A6548,6583 have strong ex- 
cess that can not be explained by the contribution 
from stellar winds and supernovae. We suggest that 
in this case the ULX is an additional source of me- 
chanical power, comparable with the total luminosi- 
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Figure 5: Integral spectrum of MF16 (NGC6946 ULX-1 nebula) taken with SCORPIO. 
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Figure 6: High- excitation ULXNe in M101 (a) and NGC5204 (b) in 16" x 16" MPFS field. To the left continual 
flux maps in MPFS fields are shown overplotted with the [O in]X4-959, 5007 flux contours. Here and below tick 
marks point towards lower intensity regions. To the right integral spectra of the [O lll]-bright regions adjacent 
to the X-ray sources are shown. Chandra X-ray positions are shown by circles. 



ties of the brightest ULXNe like MF16. The residual 
emission line luminosities are given in brackets in ta- 
ble n 

M51 X-7 is located at the outskirts of a bright star 
cluster (figure [9]), clearly identified with a young mas- 
sive cluster n5 194-839 from the list of Larsen (2000). 
It is a Super Star Cluster (SSC), or a young mas- 
sive cluster, with the integral My — — 11 I ?09 and 
U — B = — 0™81, indicating a rather young object 
of T ~ 12 Myr. However, our spectral fitting with 
StarBurst99 models gives much larger cluster age 
T ~ 60 ± 15 Myr and A v ~ m 5 ± 0.1. Most likely, 
stellar populations of different ages are seen in the 
vicinity of the X-ray source. Ha emission accompa- 
nied with bright [N n]A6583,6548 and [S n]A6717,6731 
emission lines is coincident with the X-ray source 
within the spatial resolution. As emission line spec- 
trum contains signatures of shock excitation we apply 



equation[3]for the estimate of the total power. A detail 
study of the M51 X-7 will be published in a separate 
paper. 

4. Discussion 

4.1. Metallicity and Oxygen Abundances 

In table [5] we list oxygen abundance and metallicity 
estimates for the host galaxies and the ULX positions 
in the galaxies according to Pilyugin et al. (2003). 
It was argued by Soria (2006) that low metallici- 
ties can be crucial in understanding the nature of 
ULXs. Many of the objects are indeed found in sub- 
solar metallicity environment, most likely at about 
20-50% of solar metallicity. However, in our study of 
the star forming region associated with NGC7331 X- 
1 we find that solar-metallicity models provide better 
fit to the cluster spectrum that subsolar metallicity 
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Figure 7: Integral spectrum of the HoII X-l nebula (Lehmann et at, 2005). Only part of the spectrum in 
4500-5500A spectral range is shown. 
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Figure 9: M51 X-7 16" x 16" field observed with the 
MPFS. Stellar continuum integrated in the total spec- 
tral range is shown with the Ha countours overplotted. 
Chandra ULX position is marked by a white circle. 

models with Z = OAZq and less. 

In HII regions it is easier to estimate oxygen and 
nitrogen abundances rather than [Fe/H]. We can es- 
timate oxygen abundance for some of our objects us- 
ing the method of Pagel (1992), based on the similar- 
ity between oxygen and hydrogen ionization poten- 
tials. The method is based on the three-level atomic 
solutions by McCall (1984), allowing reliable calibra- 
tions of the ion abundances from their forbidden lines 
relative strengths. Oxygen abundance estimates use 
two relations for two ionization states of the element: 
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+ 6.174 



lg 



'[OJJ]A3726 + A3720 
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[OJIJ]A4959 + A5007 
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1/2' 



+ 



e<-2 



')■ 



where 
and ti 



n e is the electron concentration in cm . t% 
are the electron temperatures (in 10 4 K units) 
of the gas emitting in [O III] and [O n] lines, corre- 
spondingly. ^3 can be obtained from the characteristic 
line ratio /[ O ///]A4959+A5007Mo///]A4363 sensitive to 
electron temperature (Osterbrock, 1974). Because of 
the close ionization potential ti must be similar to the 
electron temperature estimate t n made from another 
characteristic line ratio, /[tv//]A6583+A654s/^[w//]A5755- 
Usually, both [O n] and [O III] lines are present in a 
single H n region. Total abundance can be found by 
adding the ion abundances: 



12 + lg — 



12 + lg — 



91 

H+ 



Q++ 
H+ 



Unfortunately we do not detect [O n]A3727 line in 
our spectra because of the limited spectral range for 
MPFS and sensitivity decrease for SCORPIO in the 
near UV range, so some additional assumptions will 
be needed. 

For MF16 we have temperature estimates ts = 
1.9 ± 0.1 and t 2 ~ t n = 1.5 ± 0.2, and n e ~ 500cm- 3 
derived from the sulfur doublet [S n]A6717, 6731 
components ratio. For MF16 oxygen abundance can 
not differ much from the solar, because both very 
bright [O in] and [O i] lines are well detected. For 
a very broad range of CLOUDY models (success- 
fully predicting most of the observed line ratios) 

%)//]A3726+A3729 /%?///] A4959+A5007 ~ 1- Using this 

estimate together with the directly observed value 
^[o///]A4959+A5007/^///3 — 9, the temperature and den- 
sity estimates one finds for MF16 12 + lg ^ ~ 8.5, i. 
e. practically solar abundance. 



4.2. UV Sources 

Point-like optical counterparts of ULXs arc usually 
found as blue stars with My ~ — 5 -f- — 8 m (Terashima 
et al., 2006). On the other hand, there are strong indi- 
cations that at least some of the ultraluminous X-ray 
sources must be also ultraluminous EUV sources (in 
the range 20A < A < lOOOA responsible for hydrogen, 
oxygen and helium ionization). In the framework of 
the two most popular models (IMBH binaries and su- 
percritical accretors) , we can expect two explanations 
for the bright EUV radiation: 

- Very massive (M > 10 4 M Q ) IMBH accreting 
at several percent of its Eddington level. Less massive 
IMBHs have too hard standard disk spectra. Note 
that the donor star contribution in the EUV region 
is most likely negligible. 

- Supercritical accretion disk (SCAD) with stel- 
lar mass black hole and UV-radiating wind photo- 
sphere, like what we observe in SS433 (Fabrika, 2004). 

IMBHs in binary systems are expected to have 
nearly-standard disks (because the accretion rates are 
most likely in the range 0.01 — 1.0) with low outer 
temperatures. If the accretion disk in such system is 
tidally truncated by the secondary, its radius will be 
some fraction (roughly 0.5 of the black hole Roche 
lobe radius, practically equal to the binary separa- 
tion, for example sec Blondin (2000)) of the binary 

1/3 

separation, that is a ~ (Mbh/M*) 1 . This leads 
to at estimate of the outer disk temperature: 
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lORr. 
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\W 3 M Q 
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In UV and optical ranges IMBH accretion disks have 
approximately power-law spectra with F v oc v 1 ^ 3 
(Shakura & Sunyaev, 1973). 

According to Poutanen et al. (2006), the temper- 
ature of the outer photosphere of a supercritical disk 
wind can be expressed as functions of the black hole 
mass in solar units m and mass accretion rate in Ed- 
dington units rh as: 

J 0.2keV vnT^l^mT 1 for v oc r -1 / 2 , , 
ph = { 0.8keV m" 1 / 4 ^ 3 / 4 for v = const, W 

where v is the radial velocity of the outflowing wind 
forming inside of the spherization radius, v oc r -1 / 2 
corresponds to the case when the velocity is every- 
where proportional to the virial, and v = const is 
the case of a relatively fast wind, m is the black hole 
mass in Solar units, rh is the mass accretion rate in 
Eddington units. 

For accretion rates rh ~ 10 — 100 maximum of 
the photosphere radiation falls in a spectral range of 
hundreds of angstroms. Expected total photosphere 
luminosity is ~ 10 39 -j- 10 40 ergs s _1 . 

In figure [10] we present spectral energy distribu- 
tions (SEDs) for multicolor accretion disk (Shakura 
& Sunyaev, 1973; Mitsuda et al., 1984) models for 
IMBHs and SCAD photospheres with different black 
body temperatures. The NGC6946 ULX-1 SED in- 
ferred from optical and X-ray observations is also pre- 
sented. The SED of NGC6946 ULX-1 looks fairly flat 
compared with blackbody and MCD models. This can 
be explained qualitatively by a supercritical accretor 
system spectrum (Fabrika et al., 2007): the disk fun- 
nel (~ 10 39 — 10 40 ergs s _1 in X-ray band), the wind 
photosphere (~ 10 39 - 10 40 ergs s -1 in EUV/UV) and 
the donor star (~ 10 38 ergs s _1 in the optical/UV). 

If ULXs are indeed supercritical accretors one may 
expect a wide range of X-ray luminosities and pho- 
tosphere temperatures, resulting in large variety of 
UV/EUV properties practically independent of the 
X-ray spectrum. For rh ~ 100 a bright EUV source 
and a photoionized nebula are likely to appear. For 
rh < 10, however, one can expect only a compact HII 
region with relatively bright high-excitation lines like 
He IIA4686 and [O m]A5007. This is probably the case 
for HoII X-l and M101 P098. Very large accretion 
rates (rh > 1000) result in a very soft photosphere 
emission combined with X-ray ionization and strong 
wind producing a large shock-powered nebula. 

As can be seen from figure [TO] ULXs are reachable 
sources for GALEX (Martin et al., 2003). For sources 
with lower absorption even UV spectroscopy will be 
possible in pointing observations. The main predicted 
difference between SCADs and IMBHs in the UV is 
in the difference of the spectral slopes in these two 
models. We also expect SCADs to be more divergent 
spectral properties in the UV range. 



5. Conclusions 

We confirm that ULXs belong to young stellar popu- 
lation. Some of them are certainly formed in star clus- 
ters and associations of 5-10 Myrs old. Most of ULXs 
are surrounded with nebulae. Studying the nebulae 
gives us a powerful tool to investigate the nature of 
ULXs. At least 4 (group A) of 8 ULXNe studied in 
this paper show high flux ratio [O in] A5007/H/3 = 
3-^7 (table [2]), that requires a strong EUV source 
(ultraluminous UV sources, UUVs) to ionize the neb- 
ula. Other 4 of the 8 ULXNe (group B) are most likely 
shock-ionized. 

There is a ULXN in the group B, the nebula 
of NGC7331 X-l. Abolmasov et al. (2007b) have 
found parameters of the young super-stellar cluster 
(SSC) hosting the ULX and calculated the nebula 
spectrum using CLOUDY photoionized by the SSC 
in NGC7331. Residual emission line spectrum (val- 
ues in brackets in table [2]) , obtained after subtrac- 
tion of the model spectrum from the observed spec- 
trum, it shows both enhanced shock-excited lines (like 
[S n]A6717,6731) and [O m]A4959,5007 doublet rela- 
tive to the total emission-line spectrum of the cluster 
HII region. The residual spectrum of the nebula in 
NGC7331 is quite similar to spectrum of MF16, but 
about twice as fainter. NGC7331 X-l residual emis- 
sion may be a fainter equivalent of MF16. 

Signatures of high excitation (He II, [Fe III], 
[O in] lines) are present in most of the spectra. Low- 
excitation shock-excited lines are unusually bright in 
the spectra of ULNXe as well. We isolate then a 
group A (HoII X-l, NGC6946 X-l, M101 Po98 and 
NGC5204 X-l) requiring a strong EUV source to ion- 
ize oxygen and a group B (HoIX X-l, IC342 X-l, 
NGC7331 X-l and M51 X-7) having higher ratios 
[SII]A6717,6731/H/3 = 2.3-4 and requiring shock 
excitation of their nebulae. We note that in fact all 
the ULXNe studied in the paper show the shock ex- 
citation, because their ratios [SII]A6717, 6731/Ha 

().:! i.i 

We may make (preliminary) conclusions: 

- All the ULXNe show signatures of shock exci- 
tation, however the local ISM environment can also 
change the appearance of a ULX nebula. 

- Half of the ULXNe (the group A) require 
additional strong EUV continua (ultraluminous UV 
sources, UUVs) to ionize the nebulae. 

However, this isolation of A and B groups may be 
done only in average. For example, IC342 X-l neb- 
ula (the group B) shows that the central parts of the 
nebulae is bright in high-excitation lines. In fact all 
the ULXNe except two (M51 X-7 and NGC7331 X-l 
residual emission) show either Hell A4686 emission 
or enhanced [O in] A4959, 5007 / H(3 ratio which in- 
dicate the presence of hard ionizing source with the 
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Figure 10: Spectral energy distribution from the X-rays to the optical range derived for NGC6946 X-l. Optical 
source (star d, the V-band flux) is indicated by an asterisk, the upward arrow is the value corrected for the 
Galactic absorption (it is a lower limit because of possible intrinsic absorption). Two diamonds with bars 
ranging towards higher energies are the estimates for the EUV luminosity based on the optical emission-line 
luminosities of MF16 (HI and Hell lines, the symbols mark HI and Hell Lyman edges). The X-ray spectrum 
is represented by the best-fit model for the Chandra data given by Roberts & Colbert (2003). In two GALEX 
bands the limiting fluxes are shown with the S/N = 3 detection limit for exposure time t = 10 5 s. (a) MCD 
models for black holes accreting at O.OlMEdd are shown with solid curves. Dashed lines show both a power-law 
interpolation between hydrogen edge and V-band, and the same spectrum, but reddened by Galactic absorption 
(Cardelli et al.(1998) reddening curves were used), (b) Black body models with L = 10 39 ergs s~ 1 are presented, 
together with the optimal blackbody model for the UV radiation obtained in CLOUDY modelling. Dotted curves 
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luminosity ~ 10 38 ergs s _1 . Most likely shock waves, 
X-ray and EUV ionization act simultaneously in all 
the ULXNe, but the relative importance of the power 
sources varies significantly. 

In some ULXNe radial velocity gradients 
(~100km/s) were directly detected (Dunne et 
al.(2000), Pakull et al.(2006), Fabrika et a!. (2006) 
and references therein). We make an additional 
conclusion: 

- ULXs must produce strong winds and/or jets 
powering their nebulae with > 10 39 ergs s -1 . This is 
consistent with the suggestion that ULXs are high- 
mass X-ray binaries with the supercritical accret 
of the SS433 type. From the other hand, accreting 
IMBHs with standard disks are unlikely to have jet 
or wind activity. 

Some of ULXs are bright ultraviolet sources with 
hydrogen- and helium-ionizing luminosities varying 
from ~ 10 38 to ~ 10 40 ergs s _1 . This property can be 
still explained by a population of IMBHs with com- 
paratively high (~ 10000M©) masses or stellar-mass 
supercritical accretors (SCAD) like SS433. Future UV 
(Galex) observations may help in solving of this task, 
we predict different spectral slopes in the UV spectral 
ranges in the IMBHs/SCADs models. 
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